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Distinguishing the signals due to scattering of WIMP dark matter off of nuclear targets from those 
due to background noise is a major challenge. The Earth's motion relative to the galactic halo 
should produce halo-dependent seasonal modulation in the event rate, but it also should produce an 
angular signal that is both far stronger and less ambiguous. Distinct patterns in the recoil spectrum 
can reflect the details of the galactic halo. We derive a new formalism to calculate angular event 
rates, and present the predicted angular signal for a variety of halo models and calculate the number 
of events needed to distinguish a dark matter signal from an isotropic background. 
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The direct experimental detection of WIMP dark mat- 
ter would have profound implications for both particle 
physics and cosmology. Unfortunately, the signature, 
excess energy deposits due to scattering on nuclei, is 
not easily distinguishable from various radioactive back- 
ground signals. A significant number of events might 
be needed before a definitive observation is believable. 
Moreover, the recoil spectrum alone is largely indepen- 
dent of many features of the galactic halo distribution, 
so that our ability to distinguish between galactic halo 
models will be limited by such observations. 

In order to distinguish and isolate the WIMP signal, 
new signatures are of great interest. It was recognized 
quite early on that the Earth's motion through the galaxy 
should induce both a seasonal variation in the overall 
event rate and an overall forward-backward asymmetry 
in any directional signal The seasonal modulation 

is problematic, however. Not only is it very small, of 
the order of (Zw^/whaio) 2 ~ 0.03 — 0.05, but as has be- 
come clear as a result of recent claimed dark matter de- 
tections, the backgrounds themselves are likely to have 
seasonal modulation, via modulations in the such things 
as cosmogenic production of radon, etc. 

The directional signature is, a priori, preferable in al- 
most every way. The forward backward asymmetry can 
be large, 0(i>©/i>haio) ~ 1, and backgrounds are unlikely 
to reproduce this signature. Nevertheless, in spite of its 
intrinsic interest and probably because detectors with an- 
gular sensitivity are not yet on line, no systematic study 
of possible angular signals has yet been performed. But 
because of the difficulty of building such detectors, it is 
worthwhile to examine in advance what might be possible 
to learn with directional sensitivity. 

Indeed, one of the virtues of this approach is also one 
of its drawbacks. Angular resolution would provide sensi- 
tivity to the detailed features of the galactic halo WIMP 
distribution. However, because we do not know the na- 
ture of this distribution in advance, we cannot assume 



its form in advance in order to search for a signal. In 
fact, one might imagine that existing uncertainties are so 
great that one might not be able to obtain any unam- 
biguous limits on the basis of any lack of asymmetry in 
the observed signals. It is clear that the estimates which 
have previously been carried out, for a spherically sym- 
metric isothermal halo, are not sufficient to determine in 
general how many events may be needed to distinguish 
signals from backgrounds. In the past decade we have 
learned that the galaxy has large asphericities. It has 
also has been argued that the halo dark matter distribu- 
tion may not resemble an isothermal gas at all || . 

In order to address these features, we have developed 
a new formalism which allows one to calculate the full 
differential event rate in detectors, as a function of both 
energy and angle, for any incident WIMP distribution, 
without assumptions of spherical symmetry or even cylin- 
drical symmetry. In addition, the full motion of the 
earth around the Sun is taken into account, not merely 
the component of the Earth's velocity tangential to the 
Sun's galactic velocity. While the perpendicular compo- 
nent is largely irrelvent for estimating angular dependent 
effects if the incident halo is spherically symmetric, this 
is not the case if this assumption is relaxed. We report 
here on our general results, leaving a detailed derivation 
to a later publication. In addition, we summarize the 
main results from our first application of this formalism. 
We have generated Monte Carlo distributions for the full 
spectrum of reported halo models, in order to perform a 
statistical analysis to determine to what extent any di- 
rectional signal might unambiguously be distinguishable 
from an isotropric background. Our results are strking. 
We find that, independent of the halo model, for energy 
thresholds and quenching appropriate to the present gen- 
eration of WIMP detectors, if fine scale angular resolu- 
tion is possible, fewer than 25 events will be needed to 
distinguish the halo signal from a uniform background 
if the signal to noise ratio is large, and fewer than 50 
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events would be required if the signal to noise ratio is 
unity. On the other hand, if only a forward-backward 
asymmetry is discernable, then between 500-3000 events 
may be required to distinguish a signal, depending upon 
the actual halo distribution. In future work we will ex- 
plore how one might distinguish between different halo 
models, and how varying angular and energy resolution 
and WIMP cross sections will effect these estimates. 

Consider a WIMP of mass m x with its direction spec- 
ified by the two angles (a, 0) incident on a fixed nucleus. 
The WIMP elastically scatters off the fixed nucleus and 
the scattering event can be described by the two angles 
(#*,£) defined in the center- of-mass frame. The final 
angle set (7, <j)) describes the direction of the recoiled nu- 
cleus in the lab. Note that only the direction of the 
recoiled nucleus is visible to the detector. Neither the 
incident nor the scattered direction of WIMP can be ob- 
served. The probability of an incident direction (a, (3) is 
given by the halo model which defines the WIMP distri- 
bution function. Once we relate the three set of angles 
discussed above using kinematics, we have only to find 
the proper Jacobian transformations in order to present 
the final event rate as a function of the visible angles 
(7,0). This is simplified by the introduction of a single 
function J(a,/3;7, <fi) defined as a scalar product of two 
unit vectors on the sphere, 

J(a, /3;7, (p) = [cos7COsa + sury sina cos(</> — /?)]. (1) 

This observation leads to the general angle dependent 
event rates dR/dfl and d 2 R/dQ dil for arbitrary distri- 
bution functions. Here Q(v,J) is the energy transferred 
to the nucleus during the collision with the WIMP and 
the is the angle along which the target nucleus recoils. 
For this work we focus on the angular distribution 

dR o- Py 
= — x 

dfl 7i< £ irm n m x 

jf °' C v 3 dvF 2 (Q) J dn a ,fif{v,a,ff)J(a,Pn,<f>)e(J) (2) 

where v m [ n is related to the threshold energy, ctq is 
the low energy WIMP-nucleus scattering cross section, 
v csc = 650 km/s is the escape velocity of the galaxy, p x 
is the local dark matter halo density, O(J) is a unit step 
function, m n is the mass of the target nucleus, and F (Q) 
a the form factor suppression for scalar interactions when 
the WIMP mass couples to the quantum numbers of the 
entire target nucleus. In this work, we take a simple expo- 
nential form for this form factor F 2 (Q) oc exp(— Q/Qo), 
where Qo — 2m 3 Ri and m n is the mass of the target 

nucleus and Rq = 0.3 + 0.91^/m„ is the radius of the 
nucleus (in femtometers when m n is in GeV). 

This formulation of the differential cross section is gen- 
eral enough to accomodate any incident WIMP distribu- 
tion function f(v, a, 0). However, when we evaluate these 



event rates, extra care is needed in the lower limit of the 
velocity integration defined by the threshold energy of the 
detector. Unlike the angle independent analysis where 
the lower limit is fixed once we have the specific threshold 
or the expected transfered energy, the limit here varies 
as the incident direction changes for a fixed detector po- 
sition. This is the main source of the variation of the 
event rate as a function of angle in the angle dependent 
spectrum. This change is also taken care of with a di- 
vision of the lower limit by the same Jacobian (g). The 
angular distributions shown have been averaged in 5 day 
bins over the Earth's motion for one year. 

We next turn to the galactic halo models. Our de- 
sire here is to span the range of realistic possibilities 
that have been explored in the literature, in order to 
get a good idea of the existing uncertainties. This might 
be an important consideration for experimentalists who 
may devote considerable time to attempting to features 
in the WIMP distribution that might not be guaranteed 
to exist. We consider the following: isothermal models, 
"Evans" models 0] which are axisymmetric and allow 
for flattening, (previously been studied in the context 
of annular modulations [El), co-rotating models with a 
net angular momentum for the halo 0, and finally the 
most exteme departure from an isothermal model, using 
phase space flows dominated by infall into the galaxy for 
non-dissipative WIMPS, the "Caustic" model f|. 

For the isothermal model we will consider three ve- 
locity dispersions, (v 2 ) — 3«q/2 where vq — 150 km/s, 
220 km/s, and 300 km/s. For the Evans model the im- 
portant parameter is the flattening, q. We consider the 
values q — 1 (cored isothermal distribution), 0.85, and 
l/\/2 (maximal flattening). For the Caustic model the 
parameters for the dark matter streams are given in table 
[1] ||. If we assume the local density of dark matter is 
0.52 GeV cm -3 , appropriate to fits to the galactic rota- 
tion curve in this model then by summing the density in 
each caustic peak we find that 60% of the local density 
is due to caustics. The other 40% we assume comes from 
an isothermal distribution with vq — 220 km/s. 

To model a rotating halo we follow a standard prescrip- 
tion @. Let f + (v,a,0) = f(v,a,0) when < a < ir/2 
and otherwise. Similarly let f—(v,a,/3) — f(v,a,/3) 
when 7r/2 < a < 7r and otherwise. The distribution for 
a rotating halo is then given by / = + + — «)/-. 
Note that (v 2 f^ = (v 2 f). The average velocity of the 

halo matches the velocity of the Sun when k = 

With these analytic distribution functions, we can uti- 
lize our general formalism to calculate the angular event 
rates, forward backward asymmetries, and annual modu- 
lations expected in detectors. A central issue mentioned 
above, and the one we focus on here, is whether exist- 
ing uncertainties in galactic halo models might invalidate 
the approach of utilizing angular resolution to extract a 
WIMP signature from an otherwise isotropic background. 
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To display event rate distributions we choose a set of 
detector parameters. We use germanium (mQ C = 73 GeV) 
as the target nucleus. Note that germanium has a 
quenching factor of 0.25, thus only 25% of the incident 
energy gets transfered to the recoiling nucleus. This is 
important for determining detector thresholds. In the fig- 
ures an assumed WIMP mass of 60 GeV was used, and 
thresholds incorporate this quenching factor. We did not 
assume a specific particle cross section, detector size nor 
efficiency. Instead we calculate the number of events re- 
quired for identification of a signal. Based on this number 
and the local halo density the detector size and efficiency 
can be determined for a given cross section. 
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FIG. 1. Angular event rate distribution for an isothermal 
model with vo = 220 km/s, and for a caustic infall model 
with parameters described in the text. The isothermal model 
falls off exponentially in the forward direction while the caus- 
tic model peaks slightly in the forward direction but is more 
isotropic due to the presence of both caustic and isothermal 
components in the model. 

We employ a maximum likelihood analysis, along with 
Monte Carlo generation of sample scattering distribu- 
tions, as follows. We define a likelihood function 

We 

£ = JJP(7i,M, (3) 

i=l 

where N e is the total number of events and P(pji,<f>i) 
is the probability of a nuclear recoil in the ji, 4>i di- 
rection based on a particular model (e.g. an isother- 
mal distribution). At the 95% confidence limit when 
log/I^/dQ — log£fl at < 1 the two distributions are in- 
distinguishable. We generate 10,000 data sets for each 
N e and demand that the log-likelihood condition is satis- 
fied less than 5% of the time. The smallest N e for which 
this occurs is the minimum number of events required to 
get a 95% detection 95% of the time. 

With the above procedure we can determine the num- 
ber of events required to distinguish the signal from a flat 



distribution, both for a pure signal and when the signal- 
to-noise ratio, S/N, is one. For S/N = 1 we replace P 
with Ptotal = APflat + (1 - X)P d R/dfi where A = 1+ g/ N - 
We can perform a similar analysis for the forward-to- 
backward asymmetry. In this case the probability func- 
tion P is given by a binomial distribution. 
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FIG. 2. Required number of signal events for isother- 
mal halos against the null hypothesis. Solid lines correspond 
to the vo = 220 km/s and dashed lines to vo = 300 km/s 
. For the former, both S/N= 1 case (upper) and no noise 
cases are shown as are both the number of signal events as a 
function of threshold and as the equivalent number of events 
required for zero threshold. For the other case only the latter 
curve, for S /N= 1 is shown. This latter curve is also displayed 
when only the forward-backward asymmetry is used to probe 
a»[) = 220 km/s isothermal halo. 

The results as a function of the threshold for an 
isothermal halo with velocity dispersions of 220 km/s 
and 300 km/s are given in Figure |^ and for a forward- 
backward asymmetry measurement for the 220 km/s 
halo. In this figure, for the 220 km/s halo case, the bot- 
tom sets of lines give the number of signal events required 
as a function of the threshold. Fewer events are required 
to distinguish a halo from a flat distribution in this case 
at higher threshold because the anisotropy of the distri- 
bution becomes more dramatic as higher scattering en- 
ergies, and thus higher incident velocities relative to the 
Earth, are required. However, while fewer events are re- 
quired as the threshold is increased, even fewer events 
are expected, for a given detector size, as the threshold 
is increased. Thus, one does not win by simply raising 
the theshold. To make this clear, the upper curves show 
the number of events (proportional to the size of the de- 
tector) that would be needed to be observed in a zero 
threshold detector in order to obtain the requisite num- 
ber of events at each non-zero threshold. This curve is a 
monotonically increasing function of threshold. Clearly, 
the lower the threshold, then, the smaller the detector, 
and the smaller the total number of scatterings required 
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in the detector in order to distinguish the distribution 
from a flat distribution. 

The upper two sets of curves for the 220 km/s halo case 
in the figure represent the same requirement, but now in 
the case when the signal to noise ratio, S/N = 1 at each 
value of the threshold. The number of events plotted are 
the signal events. The noise events would lead to twice 
as many total events actually being detected. For the 
300 km/s case, the forward-backward asymmetry, and 
the the co-rotating halo, Caustic, and Evans models (fig. 
3) only the latter S/N = 1, zero threshold-equivalent 
number of events are shown. 
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FIG. 3. Required number of signal events for the Caustic 
model, a co-rotating, and an Evans model (with q = \/\/~2) 
against the null hypothesis, as in figure 2. Only S/N = 1, zero 
threshold-equivalent number of events are shown for these 
models. 

Perhaps the most suprising result of our analysis is 
that, independent of the existing uncertainties in the 
halo distribution, that less than 50 detected events for 
a signal to noise ratio of unity would be required in or- 
der to unambiguously distinguish a halo induced signal 
from that due to a flat background, assuming fine grained 
angular resolution. If, however, only forward/backward 
sensitivity were available, the requisite number of events 
increases dramatically for non-isothermal distributions. 
This number rises to 500 for the co-rotating case, and 
3000 for the Caustic case. Next, we note the suprising 
result that the Caustic distribution, which in principle is 
the most distinctive halo distribution of all, requires the 
greatest number of events in order to be distinguished 
from a flat distribution. There are two reason for this. 
First, because this distribution is dominated by infall 
components, the induced anisotropy in the pure cautic 
distribution is in the opposite direction to that predicted 
in a standard isothermal model. Adding a isothermal 
component then tends to cancel the anisotropy induced 
by the Caustic flows. In addition, because there is no ex- 
ponential tail for the caustic, the aniostropy tends to rise 



linearly in the forward instead of exponentially (in the 
backward direction) as in the pure isothermal models. 

For the same reason it is clear that the Caustic distri- 
bution should be more easily distinguishable from these 
other distributions than it is from a flat background. 
The formalism we have presented here will allow quan- 
titative estimates for any observables to be derived for 
any incident halo distribution. Indeed, in future work 
we plan to investigate how many events will be required 
in detectors in order to distinguish between the various 
halo models, as well as incorporating additional detec- 
tor features such angular resolution and energy detec- 
tion uncertainties as well as the theoretical distinctions 
between spin-dependent and spin-independent scattering 
cross sections. Our initial results, however are encour- 
aging. Independent of halo uncertainties, good angular 
resolution would allow an anisotropic halo signal to be 
differentiated from a flat background signal with a frac- 
tion of the number of events that would be needed to 
probe for annual modulation, or other expected spec- 
tral features due to a WIMP halo in the galaxy. If a 
forward-backward asymmetry is all that can be detected, 
the number of events required is model dependent, and 
can be quite large, however. 

This work was supported by the DOE. 
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TABLE I. 


Velocity flows of dark matter from the Caustic 


model (for h 


= 0.75, 


£ = 0.28, jmax 


= 0.25 Sikivie model). 


Velocities are 


given in 


the rest frame of the galaxy. 




Flow 


P a 


v x b 




v z b 


1 


0.4 


140 





±600 


2 


0.9 


250 





±500 


3 


2.0 


350 





±395 


4 


6.1 


440 





±240 


5 


9.6 


440 


±190 





6 


3.0 


355 


±290 





7 


1.9 


295 


±330 





8 


1.4 


250 


±350 





9 


1.0 


215 


±355 





10 


1.1 


190 


±355 





a In units of 10~ 26 g cm' 3 



b velocities shown are in km s 
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